Metal-terminated Graphene Nanoribbons 
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We have investigated structure, electronic, and magnetic properties of metal-terminated zigzag graphene 
nanoribbons (M-ZGNRs) by first-principles calculations. Two families of metal terminations are studied: 1) 
3d-transition metals (TM) Fe, Co, and Ni; and 2) noble metals (NM) Cu, Ag, and Au. All systems have 
spin-polarized edge states with antiferromagnetic (AFM) ordering between two edges, except Co-ZGNRs and 
Ni-ZGNRs which exhibit negligibly small energy differences between AFM and ferromagnetic (FM) states with 
the given ribbon width. In the AFM state the TM terminations transform semiconducting ZGNRs into metallic 
ones, while the band gap remains in ZGNR with NM terminations. Ferromagnetic states of M-ZGNRs with TM 
terminations show a high degree of spin polarization at the Fermi energy. We predict a large magnetoresistance 
in Fe-ZGNR junctions with a low, uniform magnetic switching field. 

PACS numbers: 73.22.Pr, 75.47.-m, 71.15.Mb 



I. INTRODUCTION 

Giant magnetoresistance (GMR)''^ and tunneling magne- 
toresistance (TMR) effects^"^ in spin-valves and magnetic tun- 
nel junctions have been the focus of spintronics studies in 
the past two decades. Conventional GMR and TMR devices, 
based on inorganic multilayered structures whose critical ele- 
ment is a sandwich of two magnetic layers separated by non- 
magnetic metallic or insulating layers, exhibit a large change 
in conductance when the orientation of the magnetic layers 
is changed from parallel (P) to antiparallel (AP). The mag- 
netic layers are usually made of transition metal ferromagnets, 
which gives a high spin polarization at the Fermi energy. 

Graphene nanoribbons (GNRs) with armchair (AGNR) or 
zigzag (ZGNR) edges**"'^, a novel organic material system 
produced by cutting graphene along two high-symmetry crys- 
tallographic directions, have recently attracted attention in 
spintronics. ZGNR is predicted to be a semiconductor with 
spin-polarized electronic states localized on the two edges 
that are individually ferromagnetically (FM) ordered but an- 
tiferromagnetically (AFM) coupled to each other through the 
graphene backbone'^. Son et al}'^ show that ZGNRs become 
half-metallic when an external transverse electric field is ap- 
plied across the edges. Kim et alM propose a spin-valve 
based on a ZGNR connected by two ferromagnetic electrodes 
using a local magnetic field to switch the magnetization of 
one of the electrodes. Very recently a magnetoresistance de- 
vice based on the transformation from the semiconducting 
state to the metallic state of the ZGNR has been proposed 
by Munoz-Rojas et al}^; an estimated 200 T magnetic field 
is required to switch the AFM to the FM state at room tem- 
perature. Although these work demonstrate convincingly that 
spin-polarization and magnetoresistance can be realized us- 
ing ZGNRs, the large magnitude of the electric field required 
to close the bandgap, the need of local magnetic control or a 
huge magnetic field, and small spin correlation length due to 
the lack of magnetic anisotropy in ZGNRs'^ are serious lim- 
iting factors for future applications. Moreover, the predicted 
magnetic states in ZGNRs have not been experimentally ob- 
served so far. 

Most of the theoretical studies on GNRs for spintron- 



ics applications so far only address ribbons with hydro- 
gen terminationsii^. Recent experiments have shown that 
graphitic surface curvature can be used as constraint and guide 
in which metal clusters aggregate to form linear islands'**^^. 
The presence of metal atoms adsorbed on the graphena2£i2^ 
or GNR^'* changes its electronic and magnetic properties and 
thus is certainly of great interest in physics. 

In this paper we show that hybrid metal-terminated ZGNRs 
(M-ZGNRs) can be excellent candidate for spintronic appli- 
cations. We studied various M-ZGNRs with terminations of 
3d ti-ansition metals (TM) Fe, Co, Ni and noble metals (NM) 
Cu, Ag, Au using first-principles calculations. Compared with 
hydrogen terminated ZGNR (H-ZGNR), the TM and NM ter- 
minations lead to drastically different and unexpected effects 
in the structural, electronic and magnetic properties of the rib- 
bon. A very large spin polarization is predicted for M-ZGNR 
with TM terminations in FM state. We propose a new type 
of junction with a low, uniform magnetic switching field that 
utilizes the electronic properties of M-ZGNRs, for which a 
magnetoresistance value as high as 200% can be obtained. 



II. COMPUTATIONAL METHODS 

The electronic structure calculations are performed using 
density functional theory (DFT) implemented in the PWSCF 
code2^. Each M-ZGNR is simulated within a supercell geom- 
etry containing two pristine ZGNR unit-cells with 40 C atoms 
and 4 metal atoms at two edges. Ultrasoft pseudo-potentials 
with kinetic energy cutoff of 500 eV are employed in all sim- 
ulations. For the exchange and coiTelation functional we used 
the Perdew-Burke-Ernzerh generalized gradient approxima- 
tion (GGA)^^. A vacuum layer of 15 A is used between two 
edges and of 15 A between two graphene planes to prevent in- 
teraction between adjacent images. Brillouin-Zone sampling 
uses a grid of 12 Monkhorst-Pack^^ fc-points along the peri- 
odic direction of the ribbon. The optimization of atomic po- 
sitions proceeds until the force on each atom is less than 0.01 
eV/A. 

The transport calculations of the proposed Fe-ZNGR mag- 
netic junction are performed by using the TRANSIESTA 
code^^, which is based on nonequilibrium Green's function 
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FIG. 1: (Color online). Four types of metal-terminated ZGNR: (a) 
Ni-ZGNR as linear A type (LA); (b) Cu-ZGNR as linear B type (LB); 
(c) Fe-ZGNR as dimerized linear type (DL); (d) Au-ZGNR as zigzag 
type (ZZ). The atomic configuration appears on the left, and on the 
right the corresponding spin density (p = p„ - pp) isosurfaces of the 
M-ZGNRs. The red (grey) and the blue (dark) isosurfaces denote 
±0.007 |e|/A^, respectively. The dashed box shows the edge view. 



formalism and DFT as implemented in the SIESTA package^. 
In order to reproduce the correct band structure of the lead, 
it is necessary to use double-^ basis with polarization for Fe 
atoms, while single-^ basis with polarization is sufficient for 
C atoms. The equivalent plane-wave cutoff for the real-space 
grid of 150 Ry is used in throughout the calculations. Current- 
voltage characteristics within small bias are calculated from 
the transmission coefficient T{E) using / = J T (£)[/(£ - 
Hi) - f(E - iUR)]dE. The magnetoresistance can therefore be 
calculated via MR - (/p - /ap)//ap- 



III. RESULTS AND DISCUSSION 

A. Structure, electronic, and magnetic properties 

To determine the lowest-energy structures of the M-ZGNR, 
first we carried out structural relaxations for four possible ter- 
mination configurations: the linear A type (LA), the linear B 
type (LB), the dimerized linear type (DL), and the zigzag type 
(ZZ), as shown in Fig. [1] Both edges of the ribbon have the 
same configuration for all systems. The calculated ground- 
state configurations as well as energetic and electronic prop- 
erties are summarized in TableU We find that for all cases the 
metal atoms bond strongly with edge carbon atoms, though 
the bonding configurations are quite different from one an- 
other, as shown in Fig.[T] For M-ZGNRs with Co and Ni ter- 



minations, the relaxed structures are LA type. For Cu-ZGNR 
the most favorable structure is LB type, which is 0.98 eV/unit- 
cell more stable than the LA type structure, in good agreement 
with the result of Wu et al?^. The Fe atoms dimerize at the 
edge of Fe-ZGNR due to the Peierls distortion^^. The most fa- 
vorable structures for Ag-ZGNR and Au-ZGNR are ZZ type, 
with a metal-metal distance larger than the length of pristine 
ZGNR unit-cell. The binding energies £/, of the metal atom 
in M-ZGNRs, defined as Eb = (£m-zgnr - £zgnr)/4 - £^°™, 
and the formation energy of the metal chain, K*""",. = 

C'-' ' tormation 

£d.ain/2 - £^°™, are shown in Table. HI The large differences 
between Eb and £fo™ation indicate a large direct binding be- 
tween the metal and carbon atoms at the ribbon edge. Bader 
analysis based on the real-space-charge density""' for the M- 
ZGNRs shows charge transfer from metal atom to C atom 
for all cases. The amounts of charge transfer for Au-ZGNR 
and Ag-ZGNR are 0.08 and 0.19 |e|/metal-atom, respectively. 
Such a difference can be important for nanocatalysis. 

For each M-ZGNR considered, the ground state has FM 
spin ordering of the metal atoms at each edge, while AFM 
order between the two edges is favored over FM, similar to 
H-ZGNRs and ZGNRs without H-passivation'^. AFM cou- 
pling of the two edges is energetically more stable due to the 
indirect RKKY exchange interaction as a result of the bipar- 
tite lattice"", since the metal terminations always stand at the 
opposite sublattices at the two edges. For Fe-ZGRNs, the en- 
ergy difference AE between the FM and AFM states increases 
with decreasing ribbon width^^^ due to increased interaction 
between enhanced spins in inner C sites inside the ribbon. 
Without the ribbon host the two Fe chains can only show di- 
rect FM exchange coupling through vacuum, and the calcu- 
lated AE of the two Fe chains becomes less than 0. 1 me V/A 
when the distance is larger than 8 A. Nonetheless, as shown 
in Table. H] the ^E for all M-ZGNRs are less than that of H- 
ZGNR, indicating a relatively small effective magnetic cou- 
pling between two edges, especially for TM terminations. 

In general, the total magnetic moment of the system comes 
mostly from the metal atoms and their nearest-neighbor C 
atoms at the M-ZGNR edges. For TM-ZGNRs, the edge C 
atom presents magnetization antiparallel to the nearby TM 
atom. For Co-ZGNR and Ni-ZGNR in particular (see Fig. 
[TJ the magnetic moments of the ribbon become more local- 
ized on the metal atoms at the edges compared to the inner 
C atoms in the ribbon, which leads to a negligibly small ex- 
change energy shown in Table. U For NM-ZGNRs, the metal 
atom shows magnetization parallel to that of the neighboring 
C atom at the edge, of order 0.01 jUB/metal-atom, mainly due 
to the proximity effect. The large difference between magnetic 
moments in TM and NM terminations reflects their different 
physical origins. In TM-ZGNRs, the relatively large moments 
of TM atoms come from unfilled tZ-shells and the TM-C co- 
valent bond reduces the number of unpaired ^/-electrons, giv- 
ing a reduced moment compared to the freestanding TM atom 
or monoatomic chain. For Cu, Ag and Au terminations with 
filled li-shells, the s valence electron saturates the dangling cr 
bonds similar to the hydrogen terminations in H-ZGNRs. The 
proximity-induced local magnetic moments for Cu, Ag, Au 
atoms are negligible, but the total edge moments of about 0.3 



TABLE I: Structural, energetic and electronic properties of M-ZGNR. Columns show: metal termination, structures, metal-metal (cIm-m) and 
metal-carbon (^m-c) bond length (in A), the binding energy Ei, and the formation energy of the metal chain Ef^^^^^.^^^ (in eV per metal atom), 
energy difference AE = EpM - £^afm between the ferromagnetic and the antiferromagnetic state (in meV/unit-cell), charge transfer of the metal 
atom (in |e|) from Bader analysis, total magnetic moment mtot of the ribbon in the FM state (in Bohr magneton //b per edge termination), local 
magnetic moment of the metal atom \iy^ and its nearest-neighbor C atom yiQ (in /jb), numbers of conduction bands crossing the Fermi level 
(wafm: the AFM state; ;ipj^ and rfpi^'- majority and minority-spin in the FM state), and spin polarization PpM of the ribbon at the Fermi energy 
in the FM state. Results for H-terminated ZGNR are also listed for comparison. 
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FIG. 2: (Color online). Energy band structures (left) and density of 
states (DOS, right, in states/eV) for the M-ZGNRs. (a): Fe-ZGNR 
in AFM state; (b): Fe-ZGNR in the FM state; (c) Au-ZGNR in the 
AFM state; (b): Au-ZGNR in the FM state. The dashed (red) and 
solid (blue) lines show the majority and minority spin, respectively. 
In the AFM state, majority and minority bands (DOS) coincide. 



ZGNRs are all semiconductor. We show the band structures 
and density of states (DOS) of Fe-ZGNR and Au-ZGNR in 
Fig.|2] For Fe-ZGNR, both AFM and FM states are metallic. 
In the FM state, more flat bands of minority spin cross the 
Fermi level compared to those of the majority spin, resulting 
in a high spin polarization at the Fermi level, as can be seen in 
Fig.|2](b). A detailed analysis of the projected DOS indicates 
that the states near the Fermi level are mainly from the Fe d 
orbital. Spin polarization PpM can be calculated as PpM = 
Wa - Ni3\/(Na + Np), where A^„ {Np) represents the DOS of 
majority (minority) spins at the Fermi energy. The f fm of 
the Fe-ZGNR at E^ is about 82% which is much higher than 
that of bulk Fe^"*'^^. This implies that in the ballistic transport 
regime, the conductance of electrons with one type of spin is 
larger than the opposite spin. Similar to Fe-ZGNRs, the FM 
states of Co and Ni terminated ribbons also show high spin 
polarization, 78% and 55% respectively. 

However, all studied NM-ZGNRs, have a bandgap compa- 
rable to that of the hydrogenated ribbon (0.32 eV) in the AFM 
ground-state. The band structure and density of states of Au- 
ZGNR as a typical case for NM-ZGNRs are displayed in Fig. 
|2](c) and (d). In the FM metallic state the Au-ZGNR shows a 
large spin polarization only at the energies slightly away from 
the Ef. We also note that the band structure of NM-ZGNR is 
very sensitive to the edge structure. The metastable LA type 
of Au-ZGNR is metallic for both FM and AFM states. 



B. Metal-terminated graphene nanorlbbon junction 



jUB/metal-atom are very close to that of H-ZGNR. Although 
Mermin- Wagner theorem.^^ excludes long-range order in ID 
spin lattice model at finite temperature which would limit the 
dimension of devices based on H-ZGNRs to several nanome- 
ters at room temperature'^, improvements can be expected in 
systems of M-ZGNRs by increasing the magnetic anisotropy 
as well as the spin correlation length. 

The ground-state TM-ZGNRs are aU metalUc and the NM- 



Our calculations of truncated Fe-ZGNRs provide an inter- 
esting result: if the transverse armchair edge is terminated by 
Fe, the strong magnetic coupling will align all edge spins in 
the same direction. Such highly spin polarized FM states can 
be useful in spintronics; in particular, a junction consisting 
of two TM-ZGNR leads separated by a spacer can operate as 
a magnetoresistance device. Fig. [3] (a) is a prototype multi- 
layer junction stacking three layers of Fe-ZGNRs. The center 
spacer can be a nonmagnetic metal such as Cr or Cu, since 




FIG. 3: (Color online). (a) Schematic figure of a multi- 
layer Fe-ZGNR/spacer/Fe-ZGNR junction prototype with antiparal- 
lel ground-state configuration, (b) Calculated ground-state (AFM) 
spin density isosurfaces of the Fe-AGNR with a width of 7.2 A. 



the behavior of the AFM exchange coupling between mag- 
netic layers through a nonmagnetic spacer is well-known^ 
and a GMR junction can thus be made; one can also use a 
small piece of non-conducting AGNR connected to the two 
leads via Fe chains, as the spacer which will provide a mag- 
netoresistance junction. The key issue is the AP configura- 
tion being the ground state through the spacer: the junction 
can be switched to the P configuration under a uniform exter- 
nal magnetic field applied to both leads, without the need for 
a local magnetic field. A conductance difference between P 
and AP is in general expected. Our calculations show that the 
Fe-terminated two armchair edges favor either AFM or FM 
configurations depending on the width of the spacer. In Fig. [3] 
(b) we show an Fe-AGNR with a width of 7.2 A which has an 
AFM ground-state that is 1 . 1 me V/A lower than the FM state. 
AFM coupling of the two armchair edge states of the spacer 
can lead to an AP ground-state magnetic coupling between the 
two Fe-ZGNR leads. This energy difference can be tuned by 
increasing the thickness of the lateral Fe chain and by using 
multiple layers as shown in Fig. [3] (a). We also note that such 
a junction can be made by two Fe-AZGNR leads separated by 
a Fe-ZGNR spacer with similar function^^. 

To illustrate the possible magnetoresistance, we have per- 
formed quantum transport calculations for a model magnetic 
tunnel junction based on two semi-infinite Fe-ZGNR leads. 
A 3.1 A gap is set between two leads to represent the Fe- 
AGNR spacer The spin-dependent cuiTent-voltage charac- 
teristics for small bias (linear response region) are calculated 
from the energy-dependent transmission coefficients in the 
ballistic tunneling regime. Large transmission differences are 
found near the Fermi energy, as shown in Fig.|4](a). By com- 
paring the current-voltage characteristics between P and AP 
configurations in Fig.|4](b), we find a magnetoresistance up to 
200% within the small bias region, which is comparable to the 
TMR ratio observed in conventional Fe/MgO/Fe sandwiched 
junctions^'^. Negative magnetoresistance is also found in the 
region from to 0.2 eV, as shown in Fig.|4](c), primary due to 
the existence of several large peaks in the transmission of AP 
configuration which are very close to the Fermi energy. 

Finally we discuss the operational conditions of our pro- 
posed magnetic junction of Fe-ZGNR/Fe-AGNR/Fe-ZGNR. 
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FIG. 4: (Color online). Transport properties for a model magnetic 
tunnel junction with two semi-infinite Fe-ZGNR leads separated by 
a 3.1 A vacuum gap. (a) Energy-dependent transmission coefficients 
and (b) current-voltage characteristics within small bias region. The 
solid (red) and dashed (blue) lines represent the results for parallel 
(P) and antiparallel (AP) magnetic configuration of the two leads, 
respectively, (c) Calculated magnetoresistance (MR) as a function of 
bias voltage. 



For a lead of infinite size any external field would switch the 
magnetization of the ZGNR lead. For a junction with finite- 
size leads, the critical switching magnetic field is defined as 
B - AE/gfiBMtot, where AE is the energy difference between 
P and AP configurations corresponding to the Zeeman energy 
splitting, Mtot is the total magnetic moment of the junction in 
the P configuration, g -2 for graphene and //b =0.058 meV/T 
A junction with two Fe-ZGNR leads of length 10 nm and a Fe- 
AGNR spacer 2.4 nm in length and 7.2 A in width produces 
a AE of about 26 meV, which is comparable to room tem- 
perature kT. The junction with ideal edge termination of Fe 
atoms has a Mtot of 415 fiB, thus resulting in a critical switch- 
ing field as low as 0.54 T, which is much less than that of the 
junction proposed by Munoz-Rojas et al}^. In experiments 
the possible formation of Fe wires of finite width at the Fe- 
ZGNR lead edges would enhance the total magnetization and 
thus further reduce the effective switching field. This junction 
shows a big advantage over the devices based on H-ZGNRs 
which mainly rely on ferromagnetic electrodes as well as lo- 
cal magnetic control'^ or a huge switching magnetic field'-. 
Current nano-size lithographic techniques can provide the ba- 
sis for fabricating and patterning of such junctions. 



IV. CONCLUSION 

In summary, we have found that TM and NM terminated 
ZGNRs have completely different structures and magnetic 
properties. The TM terminations transform semiconducting 
ZGNRs to metallic, while NM termination does not affect the 
semiconducting nature of the ZGNR. ZGNRs with TM ter- 
minations and ferromagnetically coupled edges states show a 
high spin polarization at the Fermi energy. We propose a mag- 
netic junction with low, uniform switching field using the Fe- 
ZGNRs as the ferromagnetic electrodes, and quantum trans- 
port calculations indicate a large magnetoresistance is possi- 
ble in such junctions. Our results outline a roadmap to ob- 
serve the magnetic order in graphene nanoribbon and hope- 
fully stimulate experimental activities in the near future. 



Note added in proof. In an extended study we have exam- 
ined the ribbon width dependence of magnetic coupling be- 
tween two edges of metal-terminated ZGNRs and AGNRs us- 
ing projector augmented wave (PAW) potentials within DFT, 
and the results will be published elsewhere. 
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